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We analyze the many-body properties of bilayer graphene (BLG) at charge neutrality, governed 
by long range interactions between electrons. Perturbation theory in a large number of flavors is 
used in which the interactions are described within a random phase approximation, taking account 
of dynamical screening effect. Crucially, the dynamically screened interaction retains some long 
range character, resulting in log^ renormalization of key quantities. We carry out the perturbative 
renormalization group calculations to one loop order, and find that BLG behaves to leading order 
as a marginal Fermi liquid. Interactions produce a log squared renormalization of the quasiparticle 
residue and the interaction vertex function, while all other quantities renormalize only logarithmi- 
cally. We solve the RG flow equations for the Green function with logarithmic accuracy, and find 
that the quasiparticle residue flows to zero under RG. At the same time, the gauge-invariant quanti- 
ties, such as the compressibility, remain flnite to log'^ order, with subleading logarithmic corrections. 
The key experimental signature of this marginal Fermi liquid behavior is a strong suppression of the 
tunneling density of states, which manifests itself as a zero bias anomaly in tunneling experiments 
in a regime where the compressibility is essentially unchanged from the non-interacting value. 
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I. INTRODUCTION 

Bilayer graphene (BLG), due to its unique electronic 
structure of a two dimensional gapless semiconductor 
with quadratic dispersion^, offers an entirely new set- 
ting for investigating many body phenomena. In sharp 
contrast to single layer graphene, the density of states in 
BLG does not vanish at charge neutrality, and thus even 
arbitrarily weak interactions can trigger phase transi- 
tions. Theory predicts instabilities to numerous strongly 
correlated gapped and gapless states in BLG^"^. These 
instabilities have been analyzed in models with un- 
screened long-range interactions^, dynamically screened 
long-range interactions^, and in models where the inter- 
actions are treated as short ranged"—. Irrespective of the 
form of the interaction, the instability develops only log- 
arithmically with the energy scale. However, dynami- 
cally screened Coulomb interactions have been shown to 
produce log^ renormalization of the self energy^ and ver- 
tex function'^. Such strong renormalization can result 
in significant departures from non-interacting behavior 
on energy scales much greater than those characteristic 
for the onset of gapped states. However, there is as yet 
no systematic treatment of the log^ divergences. In this 
paper, we provide a systematic treatment of the effects 
of dynamically screened Coulomb interactions, focusing 
on the renormalisation of the Green function, and using 
the framework of the perturbative renormalization group 
(RG). 

We analyze the RG fiow perturbatively in the number 
of flavors, given by = 4 in BLG. We use perturbation 
theory developed about the non-interacting fixed point, 
and calculate the renormalization of the fermion Green 
function and of the Coulomb interactions. We demon- 
strate that the quasiparticle residue and the Coulomb 
vertex function undergo log^ renormalization, while all 



other quantities renormalize only logarithmically. The 
quasiparticle residue and the Coulomb vertex function, 
moreover, are not independent, but are related by a Ward 
identity which stems from gauge invariance symmetry. 
Therefore, at log^ order, BLG behaves as a marginal 
Fermi liquid. 

We solve the RG flow equations with logarithmic accu- 
racy, finding that the quasiparticle residue flows to zero 
under RG. This behavior manifests itself in a zero bias 
anomaly in the tunneling density of states (TDOS). We 
conclude by extracting the sub-leading (single log) renor- 
malization of the electron mass, as a correction to the log 
square RG. This calculation allows us to predict the in- 
teraction renormalization of the electronic compressibil- 
ity in BLG, a quantity which is interesting both because 
it is directly experimentally measurable, and because it 
allows us to contrast the slow single log renormalization 
of the compressibility with the fast log^ renormalization 
of the TDOS. 

The structure of the perturbative RG for BLG has 
strong similarities to the perturbative RG treatment of 
the one dimensional Luttinger liquids^iii"— . We recall 
that in the Luttinger liquids, the Green function acquires 
an anomalous scaling dimension, which manifests itself in 
a power law behaviour of a quasiparticle residue that van- 
ishes on shell. In addition, the electronic compressibility 
in the Luttinger liquids remains finite even as the quasi- 
particle residue fiows to zero. Finally, in the Luttinger 
liquids, there are logarithmic divergences in Feynman di- 
agrams describing scattering in the particle-particle and 
particle hole channels, corresponding to mean field insta- 
bilities to both Cooper pairing and charge density wave 
ordering. However, when both instabilities are taken into 
account simultaneously within the framework of the RG, 
they cancel each other out, so that there is no instability 
to any long range ordered phase at low energiesii. 
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Exactly the same behavior follows from our RG anal- 
ysis of BLG, including the cancellation of the vertices 
responsible for the pairing and charge density ordering. 
However, the diagrams in this instance are log^ diver- 
gent, and even after the leading log^ divergences are can- 
celed out, there remains a subleading single log insta- 
bility. Nevertheless, this single log instability manifests 
itself on much lower energy scales than the log^ RG flow. 
Therefore, over a large range of energies, bilayer graphene 
can be viewed as a two dimensional analogue of the one 
dimensional Luttinger liquids. 

Our treatment of the log^ renormalization in BLG is 
somewhat reminiscent of the situation arising in two- 
dimensional disordered metals^^. In the latter, the log^ 
divergences of the Green function and of the vertex func- 
tion stem from the properties of dynamically screened 
Coulomb interactions, which exhibit "unscreening" for 
the transferred frequencies and momenta such that lu j 
is large compared to the diffusion coefficient. Further- 
more, the divergent corrections to the Fermi-liquid pa- 
rameters, as well as conductivity, compressibility and 
other two-particle quantities in these systems, are only 
logarithmic. This allows to describe the RG flow of the 
Green function due to the log^ divergences by a single 
RG equationi^ of the form 



(1) 



where g is the dimensionless conductance. The suppres- 
sion of the quasiparticle residue, described by this equa- 
tion, manifests itself in a zero-bias anomaly in the tun- 
neling density of states, readily observable by transport 
measurements. 



II. DYNAMICALLY SCREENED 
INTERACTION 

We begin by reviewing some basic facts about BLG. 
BLG consists of two AB stacked graphene sheets (Bernal 
stacking) . The low-energy Hamiltonian can be described 
in a 'two band' approximation, neglecting the higher 
bands that are separated from the Dirac point by an en- 
ergy gap W ~ 0.4 eV— . There is four- fold spin/valley 
degeneracy. The wavefunction of the low energy elec- 
tron states resides on the A sublattice of one layer and 
the B sublattice of the other layer. The non-interacting 
spectrum consists of quadratically dispersing quasiparti- 
cle bands E± = -iip^ /2nn with band mass m w 0.054me. 
We work throughout at charge neutrality, when the Fermi 
surface consists of Fermi points. The discrete point-like 
nature of the Fermi surface is responsible for most of the 
similarities to the Luttinger liquids. 

Although the canonical Hamiltonian has opposite chi- 
rality in the two valleys, a suitable unitary transfor- 
mation on the spin-valley-sublatttice space brings the 
Hamiltonian to a form where there are four flavors of 
fermions, each governed by the same 2x2 quadratic 



Dirac-type Hamiltonian^. We introduce the Pauli ma- 
trices that act on the sublattice space t^, and define 
T± — Ti± iT2, and p± — ipy, and hence write^^ 



H = Ho + — 



yr-^ n(x)n(x') 
2^^ |x-x'| ' 



P,(T ^ 



2m 



(2) 
(3) 



Here cr = 1, 2, 3, 4 is a flavour index, n(x) ~ '^^n^{'x.) 
is the electron density, summed over spins, valleys and 
sublattices, while the dielectric constant k incorporates 
the effect of polarization of the substrate. Note that the 
single-particle Hamiltonian Hq takes the same form for 
each of the four fermion flavors, and is thus SU{A) in- 
variant under unitary rotations in the flavor space. 

The Goulomb interaction sets a characteristic length 
scale and a characteristic energy scale ( "Bohr radius and 

r- 



Rydberg energy" ^— ■ 



10k A, 



Eq = 

KQq 



1.47 



■eV. 



(4) 



In Eq.Q, we have approximated by assuming that the 
interlayer and intra-layer interaction are equal. This ap- 
proximation may be justified by noting that the inter- 
layer spacing d w 3 A is much less than the characteris- 
tic lengthscale ao, Eq.(j4]). Within this approximation, 
the Hamiltonian ^ is invariant under SU{4) flavour 
rotations^. 

We note that for k '-^ 1 the energy Eq value is com- 
parable to the energy gap parameter W ^ 0.4 eV of 
the higher BLG bands (see Ref.Q for a discussion of 
four band model of BLG). This suggests that there is 
some interaction induced mixing with the higher bands 
of BLG. However, since a four band analysis is exceed- 
ingly tedious, here we focus on the weak coupling limit 
Eq W, where the two band approximation, Eq.®, is 
rigorously accurate. We perform all our calculations in 
this weak coupling regime, and then extrapolate the re- 
sult to Eq w 1.47 cVk"^. Since the low energy properties 
should be independent of the higher bands, we believe 
this approximation correctly captures, at least qualita- 
tively, the essential physics in BLG. Meanwhile, since W 
is the maximum energy scale up to which the two band 
Hamiltonian, Eq.®, is valid, we use W as the initial UV 
cutoff for our RG analysis. 

We wish to obtain a RG flow for the problem 1^ by sys- 
tematically integrating out the high energy modes. How- 
ever, the implementation of this strategy is complicated 



^^We have performed a unitary transformation on the Hamilto- 
nian, as outUned in Ref. 0|, to clearly manifest the symmetries. 
As a consequence, our 'valley' and 'sublattice' variables are not 
the physical valley and sublattice variables, but are linear com- 
binations thereof. 
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by the long range nature of the unscreened Coulomb in- 
teraction. Within perturbation theory, the long range 
interaction gives contributions which are relevant at tree 
level, making it difficult to come up with a meaningful 
perturbative RG scheme. Therefore, it is technically con- 
venient to perform a two-step calculation, where we first 
take into account screening within the random-phase ap- 
proximation (RPA), and then carry out an RG calcula- 
tion with the RPA screened effective interaction. We em- 
phasize that it is necessary to consider the full dynamic 
RPA screening of the Coulomb interaction, since a static 
screening approximation does not capture the effects we 
discuss below. 

The dynamically screened interaction may be calcu- 
lated by summing over the RPA series of bubble dia- 
grams, to obtain a screened interaction. The RPA ap- 
proach to screening may be justified by invoking the large 
number TV = 4 of fermion species in BLG. The screened 
interaction takes the form 



functional field integral. We have 



C/(c.,q) 



27re^ 



K|q| - 27re2n(a;,q)' 



(5) 



Here n(a;, q) is the non-interacting polarization function, 
which can be evaluated analytically^^. Here we will 
need an expression for H(a;,q) in terms of Matsubara 
frequencies w, derived in Ref. 's^, where it was shown 
that the quantity Tl{uj, q) depends on a single parame- 
ter 2ma;/g2, and is well described by the approximate 
form 



H(c^,q) = - 



Nm 



2 m 



41n^4 



2tt 



q. 

, 2m 



2 \2 



(6) 



I?^t^^,exp(-5o[^^V] ' Si[y^\^]) , (7) 



uu>^ 



-5*1 = ^^^r2C/(w,q)n„,qri_„._q + S'2. (9) 

Here the "0 fields are Grassman valued (fermionic) fields 
with fiavour (spin- valley) index a, while is a fermionic 
Matsubara frequency, F is a vertex renormalization pa- 
rameter, Z is the quasiparticle residue, and n^,q is the 
Fourier transform of the electron density, summed over 
spins, valleys and sublattices. The effective interaction 
U{uj,q) is given by Eq.®. The term S2 is included ten- 
tatively to represent more complicated interactions that 
may be generated under RG. In the bare theory, F = 1, 
Z = 1 and 52 = 0. The theory is defined with the initial 
UV cutoff Aq. Since the two band model, EqlSl is only 
justified on energy scales less than the gap W « OAeV 
to the higher bands in BLG, we conservatively identify 
Aq = W. Our main results will be independent of Aq. 

As we shall see, the RG flow will inherit the sym- 
metries of the Hamiltonian, Eq.(I2|), strongly constrain- 
ing the possible terms 82- The relevant symmetries are 
particle-hole symmetry, time reversal symmetry, SU(4) 
flavour symmetry^, and the symmetry of the Hamilto- 
nian under the transformation e**'^^i?(6'/2), where R{0) 
generates spatial rotations, R{9)p± = e^**p±. 

We will employ an RG scheme which treats frequency 
bj on the same footing as /2m^ in order to preserve 
the form of the free action Eq.® under RG. Thus, we 
integrate out the shell of highest energy fermion modes 



where A'^ = 4 is the number of fermion species. The 
dependence © reproduces H(w,q) exactly in the limits 
Lu ^ /2m and lo ^ /2ra, and interpolates accu- 
rately in between. We discover upon substituting Eq.® 
in Eq.([S]) that the dynamically screened interaction is re- 
tarded in time, but crucially is only marginal at tree level. 
It therefore becomes possible to develop the RG analysis 
perturbatively in weak coupling strength, by taking the 
limit of A^ > 1. 

Since the quantity H(cli, q) vanishes when g ^ 0, the 
RPA screened interaction ([5]) retains some long range 
character, exhibiting "unscreening" for w ^ /2m. 
This will lead to divergences in Feynman diagrams of 
a log^ character. 



III. SETTING UP THE RG 

To calculate the RG flow of the Hamiltonian, Eq.([2]), in 
the weak coupling regime, we begin by writing the zero- 
temperature partition function $ as an imaginary-time 



A' < 



f ) <A, 

2m 



(10) 



and subsequently rescale w — t- a; (A/ A'), p — )■ p(A/A')^/^, 
where z is the dynamical critical exponent^^, which takes 
value z = 2 at tree level. Because the value z = 2 is not 
protected by any symmetry, it may acquire renormaliza- 
tion corrections. However, it will follow from our analysis 
that the quasiparticle spectrum does not renormalize at 
leading log^ order, so that the exponent z does not flow 
at leading order. We therefore use z = 2 for the rest 
of the paper, which corresponds to scaling dimensions 
[uj] — 1 and [p^] = 1. Under such an RG transformation, 
the Lagrangian density in momentum space has scaling 
dimension [£] =2, and we have tree level scaling dimen- 
sions [■(/'] = 1/2 and [F] = [Z] = respectively. 

Given these tree level scaling dimension values, it can 
be seen that all potentially relevant terms arising as part 
of ^2 must involve four fermion fields. Indeed, any term 
involving more than four if) fields will be irrelevant at 
tree level under RG, and may be neglected. The terms 
with odd numbers of "0 fields are forbidden by charge 



conservation, while the quadratic terms Aijif^-ipj cannot 
be generated under perturbative RG, since they break 
the symmetries of the Hamiltonian hsted above^'^. Thus, 
the only potentially relevant terms that could arise under 
perturbative RG take the form of a four point interaction 
which may be written as 



S2 = ^J d3xd3a:'T^'^fc>i,.(a:)V'.,,(x)^I,,fc(x')V.'.;(a;'), 

(11) 

where x = (r,i), x' = {r',t'), Here T is an effective four 
particle vertex, which is marginal at tree level, the indices 
(T, a' refer to the flavour (spin- valley) of the interacting 
particles, and i,j,k^l are sublattice indices. 

The symmetries of the Hamiltonian, Eq.([21), impose 
strong constraints on the spin-valley-sublattice structure 
of the four point vertex T. Since the Coulomb interac- 
tion does not change fermion flavour (spin or valley), and 
the electron Green function is diagonal in flavour space, 
the vertex T cannot change fermion flavour. Moreover, 
the SU (4) flavour symmetry of the Hamiltonian implies 
that T does not depend on the flavour index of the in- 
teracting particles, and we may therefore drop the in- 
dices (7, cr' in Eq. PT|) . Finally, the bare Hamiltonian ([2]) 
is invariant under combined pseudospin/spatial rotations 
through e^^'^^ R{9/2). This symmetry further restricts the 
form of four point vertices in Eq. ljlll) to have sublattice 
structure Tnjj or Tijji only That is, the allowed 
scattering processes are restricted to (AA) — (AA), 
{AB) (AB) and (AB) (BA). We note that the 
processes (AB) -J> (AB) and (AB) (BA) are distinct, 
since the particles have flavour, and the interaction ([5]) 
is not short range. 

Below we obtain the RG flow for bilayer graphene, 
working in the manner of Ref . [12,] . We consider the par- 
tition function, Eq.Q, where the interaction is given by 
Eq.©. Starting from this action, supplied with ultra- 
violet (UV) cutoff Ao, we systematically integrate out 
the shell of highest energy fermion modes, Ea. dTU]) . We 
perform the integrals perturbatively in the interaction, 
Eq.(l5]). This corresponds to a perturbation theory in 
small r'^Z'^/N. We carry out our calculations to one 
loop order, and examine the renormalization, in turn, of 
the electron Green function (^Sec lIV|) . the vertex function 
r (SecEl) and the four point vertex T (Sec|Vl|. 



^■^The symmetry of the Hamiltonian may be spontaneously broken. 
However, the energy scale for spontaneous symmetry breaking is 
set by the subleading single log flowsS. and is lower than the 
energy scale for the phenomena discussed in this paper. 

•^^In that, we ignore vertices of the form T j^j^j^gd^, aabaQ^^ 
and other similar terms, which are allowed by symmetries, but 
are irrelevant in the RG sense. 



a) b) 

FIG. 1: Diagrammatic representation of self energy (a) and 
vertex correction (b) [Eas. (|14 |). (|28|) ]. Straight lines with ar- 
rows represent fermion propagator, Ea. (|13|l . wavy lines rep- 
resent dynamically screened long range interaction, Eq.©. 



IV. SELF-CONSISTENT RENORMALIZATION 
OF THE ELECTRON GREEN FUNCTION 

At first order in the interaction, the fermion Green 
function acquires a self energy S, represented dia- 
grammatically (to leading order in the interaction) by 
Figllja). A self-consistent expression for the change in 
the fermion propagator G is 

5G(u:, q) = Go(w, q)S(c^, q)Go(w, q), (12) 



6*0(0;, q) 



S(^,q) = 



iu - i/o(q)' 



(13) 



dedPp 
"(2^ 



rg[/e,pGo(£ + c^,p + q), (14) 



where S is a 2 x 2 matrix in sublattice space. 

A number of general properties of the self energy can 
be established based on symmetry considerations. It fol- 
lows from Ea. ((Ti)) that I](0,0) vanishes, since the part 
of G(e, p) which is invariant under rotations of p is an 
odd function of frequency e. Likewise, the expressions for 
diagonal entries l]/i/i(0,q) and £53(0, q), which involve 
an integral of an odd function of e, vanish on integra- 
tion over e. For the same reason, the expressions for off 
diagonal entries Tiab{^,^) and I]ba(w,0) vanish upon 
integrating the momentum p over angles. Hence, non- 
vanishing contributions arise at lowest order when the 
right hand side of Eq. ([T4l) is expanded to leading order 
in small lo and q. We obtain 

Y,Aa{^,<\) = — ]-0{uj ,ujq ,q ),(15) 

S^B(.^,q) = ^ 5(^T7^ + 0(- ,q ), (16) 

where T,aa ^^bb and TtAB = ^*ba symmetry. 

It was shown in Ref.H that idT^AA/du) and 
d'T^AB/d{q\/2m) are both log divergent, and are equal 
to leading order (see below and Sec lVIIII for alternative 
derivation). Thus the self energy can be written, with 
accuracy, as 



log^ 



E(w,q) = -jZo — Gq- 



(w,q) + In 



A 



(17) 



Here, it is understood that non- vanishing dYi/duj is due 
to the modes that have been integrated out, Eg. dTUl) . 
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Within the leading log approximation, the electron 
Green function, Ea. (jl3p . retains its non-interacting form, 
whereby the self energy, upon substitution into Eg. ([T2l) . 
can be absorbed entirely into a redefinition of the quasi- 
particle residue, as 



<5G(w, q) 



1 



sz, 



sz 



(18) 



We emphasize that the lack of renormalization of the 
mass only holds at log^ order. The subleading single log 
renormalization of the mass will be analyzed in Sec I Villi 
The renormalization of the quasiparticle residue, 
Eq. dTSl) . can be evaluated explicitly by calculating 
idTi/duj. Taking E from Eq. p^ . we write 



.9S 

duj 



j=0 



(19) 

We express the momenta in polar coordinates px = 
p cos Of, Py = psina, and straightaway integrate over 
— TT < a < TT. We further change to pseudopolar co- 
ordinates in the frequency-momentum space, e — r cos 0, 
p^ /2m = rsinS, with the "polar angle" Q < 9 < tt. Using 
the Rydberg energy i?o, Eq.Q, as units for r, we have 



duj 



dr de (sin^ 6 - cos^ 6i)r§Zo 
27r V2rsin6l- ^n(6i) ' 



(20) 



where 11(0) is the dimensionless polarization function, 
given by Eq.(l6]) with quasiparticle mass m suppressed 
and 2me/p^ — cot 6*. We note that 11(6') goes to zero 
when 6* — > 0, TT, and these zeros of the polarization func- 
tion dominate the integral and lead to the log^ diver- 
gence. Since 11(0) is even about 9 = 7r/2, the log^ con- 
tribution can be evaluated by replacing 11(6') in Eg. ([20)1 
by its asymptotic 9 ir form. 



Nm 



tan0. 



(21) 



In the region <C tt, we may approximate sin0 w 9, 
tan 6* ~ 9 and cos 6* w 1. Including a factor of 2 for the 
region ss tt, which gives a contribution identical to that 
of the region 6* ~ 0, we can express the integral Ea. (PD|) 
with logarithmic accuracy as 



i— - 2 
duj 



dr 
r 



Tr/2 



d9 



^qZq 



2t: + ^9' 



(22) 



Performing the integral over 
yields 



and assuming r ^ N'^ 



.as 2TIZq 



— In ■ 



(23) 



Integrating over A' < r < A (see Eq. (fTO|) ). we obtain 

i— = — n — ^ n- n^ -:- ■ 24) 

doj Nn^ \ 8A' A' 2 A' / ^ ^ 



We now consider an infinitesimal RG transformation. 
Defining an RG time 



Ao 
A 



S( = In 



A' 



(25) 



we rewrite the recursion relation, Eq. 



as 



*aJ = lv;^(^ + ^)^^' " = ^"^A^- 

The constant term c describes corrections subleading in 
log^, and thus may seem to be irrelevant. However, we 
shall retain it in the RG equation since it will determine 
the form of renormalization near the UV cutoff (see dis- 
cussion of TDOS in ScclVlIl). 

In our derivation of Eq. ((26l) it was assumed that our 
initial UV cutoflt Aq < N^n'^Eo/8. Such choice of Aq is 
certainly justified when N is large, which is the limit we 
worked in thus far. Better still, the condition remains 
entirely reasonable for the physical value = 4, leading 
to N'^-k^Eq/S = 24eFK~^, which is much bigger than the 
bandwidth for BLG. 

Substituting Eq. ipSl) into Ea. ([T^ . we obtain a differ- 
ential equation for the flow of the quasiparticle residue. 



az 



(C + c). 



(27) 



This equation encapsulates a one loop RG flow for the 
residue Z, describing its renormalization within a log^ 
accuracy. 



V. SELF-CONSISTENT RENORMALIZATION 
OF THE VERTEX FUNCTION T 

The screened Coulomb interaction renormalizes 
through the vertex correction, pictured in Fig[T{b). The 
RPA bubble diagrams, which have already been taken 
into account in moving from an unscreened to a screened 
interaction, Eq. ([5]) , do not contribute to renormalization. 
It may be verifled by an explicit calculation that the ver- 
tex correction in Figlljb) is given by 



ST = - 



ded^p 



(Pi 

V2ni 



(28) 



This is the same expression as for the residue renormal- 
ization [Eqs. (IT8l) .([T9 | ]. with r replacing Z, and a sign 
change. Hence, we obtain 



ar 2r3(^)z2(e) 



ae 



(e + c) 



(29) 



which is identical to the flow equation for Z, albeit with a 
reversed sign. Therefore, the product TZ does not renor- 
malize at log square order, and we can write. 



(30) 



6 



Ei-co.k,-q E,+o).k,+q E,-oj,k,-q E2+o>,k,+q 

s-0).p-q 



a 



E,.k, 



E,.k, 



b) 




FIG. 2: The renormalization of the four point vertex T 
proceeds through repeated scattering in the particle parti- 
cle channel (a) and in the particle hole channel (b), known 
as the the BCS loop and the ZS' loop in the Luttinger liq- 
uid literature^. The RPA bubble diagrams (ZS loop in the 
language of Ref . ) . which arise in the same order of per- 
turbation theory, have already been taken into account in the 
screened interaction, Eq.(l5]). 



This result is not a coincidence, since the residue Z 
and the vertex function F are not independent quanti- 
ties. The Hamiltonian, Eq.Q, is invariant under a gauge 
transformation of electron vifavefunction -0' — ipe^^, ac- 
companied by energy and momentum shifts s' — e — dtX^ 
p' = p -f Vx. This gauge invariance symmetry can be 
shown to lead to Ea. (P(I| through a Ward identity that 



relates the self-energy to the vertex function 



entity ' 



VI. RENORMALIZATION OF THE FOUR 
POINT VERTEX T 

The four point vertex T, introduced in Eg. ljlip . renor- 
malizes through the diagrams presented in Figl51Ja,b), 
which represent the repeated scattering of two particles in 
the electron-electron and electron-hole channels respec- 
tively. We follow the naming conventions used in Ref. 
in the context of the Luttinger liquid, and name these 
two diagrams, the BCS loop and the ZS' loop, pictured 
in Figllja) and Figl2ljb), respectively. In the one di- 
mensional Luttinger liquids, the two processes famously 
cancellA, so that the four point vertex does not renor- 
malize. In higher dimensions, such a cancellation is rare. 
However, the discrete nature of the Fermi surface in BLG 
results in a Luttinger liquid like cancellation of the pro- 
cesses Figl51Ja,b), as will be discussed below. 

We argued in Sec lIIII that the RG-relevant scatter- 
ing processes allowed by symmetry must have sublat- 
tice structure [A, A) -> [A, A), {A,B) {A,B) or 
{A, B) — >■ {B,A). To see the mathematical origin of such 
selection, it is instructive to explicitly write out the form 
of the electron Green function. We have 



Gaa(£,p) = 



-Zie 



£2 + (p2/2m)' 



GsB(e,p), (31) 



G^s(£,p) = — ^|^J^=G|,^(e,p). (32) 
e + (p^/2to) 



When the diagrams Figl2l[a,b) are evaluated in any chan- 
nel other than these three channels, they vanish upon 
integration over inner momentum variables, due to the 
chiral structure of the sublattice changing Green func- 
tions, Eq.(l32]). 

Similar reasoning leads to a conclusion that the 
{A, B) — )■ (B, A) vertex cannot exhibit a log^ divergence. 
As we saw above, the log^ divergences arise because the 
effective interaction f/^^p has a pole at p = and fi- 
nite e. However, the sublattice changing Green func- 
tions. Eg. ([5^ . have zeros at small p, which cancel the 
contribution of the pole in the interaction. Thus, the 
diagrams in Figl5]can only be log^ divergent if all inter- 
nal Green functions are sublattice preserving, given by 
Ea. (l5T|) . Since the process [AB) — ^ {BA) involves two 
sublattice changing Green functions, it follows that the 
integrals associated with this processes cannot be log^ 
divergent, and hence this process does not contribute at 
leading log^ order. 

Thus, at leading order, we need to consider only the 
processes {AA) -)> {AA) and {AB) {AB). More- 
over, since the interaction ([S]) does not distinguish be- 
tween sublattices, the ZS' and BCS contributions from 
FigI2ta,b) in these channels are the same. Therefore, to 
demonstrate that T does not renormalize at leading or- 
der, it is sufficient to demonstrate that there are no log^ 
divergences in the {AA) — >■ {AA) channel. 

In evaluating the ZS' and BCS diagrams (Figl2]), it will 
prove important to keep track of external momenta. The 
vertex T(i?i, i?2, i^, ki, k2, q) then represents the ampli- 
tude for the scattering process 

'0<T,A,£;i,kiV'<T',A,£;2,k2 — ^ V'cr,A,£;i+w,ki+qV'CT',/l,£;2-w,k2-q- 

Translating the ZS' and BCS diagrams in FigH into in- 
tegrals, we find the contributions 



"ZS' 
- AAAA 



^ AAAA 



J^U,^pUe-^.p-ciGAA{Ei +eM+ P) 



(27r)3 

xGAA{E2+e-oj,'k2 + p-q), (33) 
j^U,^pU,-^^p-qGAA{Ei + e, ki + p) 

(34) 



{2n) 

xGaa{E2 - £,k2 - p) 



Here, the interaction U{e,p) is defined by Eq.®, the 
Green functions are defined by Eg. ipTj) . and the integral 
goes over the shell defined by Eg. (|TOl) . 

As always in a RG analysis, we assume that the exter- 
nal frequencies and momenta are small compared to the 
internal frequencies and momenta: 



2 /2 

max I TT^-,;^ ) < A' < 

Zm Zm 



In such a case, the standard approach to handling the 
integrals over e and p involves setting the external fre- 
quency and momenta to zero at first, and restoring their 
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finite values later to regulate the infrared (IR) diver- 
gences. However, a straightforward application of this 
recipe to the integrals in Eqs. (|33| .(p4 |) proves impossible, 
because these integrals are power law divergent when all 
external momenta are set to zero. The divergence arises 
from the region near p « (which lies within the shell 
defined by Eq. ljlOp ). where the interaction is nearly un- 
screened. In this region, we have 



e,P e-a.,p-q (|p| + a|pP) (|p - q| + a|p - q|2) ' 

(36) 

with a — Ne^ /2nK. At finite q, the poles in this expres- 
sion are split apart, and thus the singular contribution of 
each pole, p = and p = q, is regularized by the inte- 
gration measure (fp so that the integrals in Eas. (|33|) . (|34)) 
remain well defined. However, when all external mo- 
menta are zero, the poles from the two interaction lines 



co-incide, and the expressions (|33|) . ([34)1 acquire a second 
order pole at p = 0. When we integrate over this second 
order pole, we pick up a power law divergence. 

Hence, if either of the ZS' or BCS diagrams existed 
in isolation, this power law divergence would indicate a 
strong (power law) instability, which would drive T into 
the strong coupling regime, where our log^ RG would 
cease to apply. However, as we will now show, the di- 
vergences in the contributions to T from the expressions 
([5^ . p4l) in fact cancel out, so that T does not flow to 
log order. To analyse the cancellation between the ZS' 
and BCS diagrams, it is convenient to add the integrands 
of Eq. ([33]) and Eq. ([34]) together before doing the integral, 
while keeping external momenta finite. Preserving finite 
external momenta ensures that the integrals Eq. (j33p and 
Eq. ([M)) are well defined. After combining the integrands, 
and denoting T^^^^ -f T^^^^ = T, we obtain 



j^U,,pU,-^^p-^GAAiEi + e, ki + p) [G, 



To simplify this expression we note that momentum q 
enters very differently in Eq. p7p as compared to other 
external frequencies and momenta Ei, E2, lo , ki, k2. 
The momentum q is needed to split the poles coming 
from the two interaction terms - if we take q to zero, the 
integral will acquire a second order pole at p = 0, leading 
to a divergence. This divergence arises from within the 
shell that we are integrating out fEg. dTUl) ). and thus the 
RG will be ill defined. In contrast, sending the frequen- 
cies and momenta Ei, E2, ui, ki, k2 to zero by applying 
Eg. ([55)1 does not cause any concern. We thus have 

T = |^C/e,pC/e,p-qG^A(e,p) 

x[GAA{e,p-q) + GAAi-e,-p)]. (38) 



{E2+e-uj,k2+p-q) + Gaa{E2 - £, ka - p)] . (37) 

I 

Interestingly, the expression in square brackets van- 
ishes identically when q = 0, since Gaa{£,p) = 
~Gyiyi(— £, — p). However, taking the limit q — > is 
potentially problematic because of the pole structure of 
C/j_pC/e_tj,p_q discussed above. Instead, we proceed with 
caution, and evaluate Eq. (|38l) at finite q, using the con- 
ditions (|35|) to simplify the analysis. 

Given what we just said, it is now easy to see why 
there is no log^ divergence in T. First, we note that the 
interaction ([5]) carries a soft UV cutoff, so the integral in 
Eq. ([38]) is UV convergent (this property of dynamically 
screened interaction in BLG is discussed e.g. in Ref.i3j). 
Hence, we can shift variables to p± — p±q/2 and rewrite 
the expression pS]) as 



de<Pp 



(27r)3^'P+ 
ded'^p 



Ue,p+U,,p_e^D{e, p+) [D{e, p_) - Die, p+)] 



(39) 



D{e,p+) + Die,p^) , Die,p+) - D{e,p^) 



[2?(e,P-)-D(£,p+)], 



r 



where we factored the Green functions as 
GAA{£,p)^iZeD{e,p), D{e,p) 



the first term in the brackets gives zero upon integration 
over p. Hence, we can rewrite the result for T, Eq.(39), 



e2 + (p2/2m)' 



We note that because T should be even under q 



(40) 
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tions for Z and F, we write the equation for Z as 



p4^2 




2 j 


' (27r)3 






2 j 


' (27r)3 



C^e.p+C/.,p_£' [^(£,P-) - ^(e,P+)] 

9 9 n 2 



(41) 



where z± = |p±p/2m. 

To extract the leading contribution at smah q, we ap- 
proximate the effective interaction as 



C/(e,p) 



n-i(e,p) 



1 



1 - 



«|p| 

27re^n(e,p) 



n(£,p) 



(42) 



From the definition of the polarization function, Eq.Q, 
we see that the approximation U ss — l/FI holds ev- 
erywhere in the shell Ea.([TO| except at p « 0, since 
n(p = 0) = 0. However, in the limit p — > 0, the ex- 
pression in brackets in Eg 1411 tends to zero because of 
the expansion z'j_ — z'^ — (p'^/m){p ■ q/2m) + 0{p*), 
which ensures validity of the approximation (|42l) . 
Hence, using Eq.®, we obtain 



T = 



p4^2 



ded?p 



AiT{Nm In 4)2 

^2 ^2 H 2 



[e-^ + zl){e-^ + zl) 



(43) 



Simple power counting shows that this integral is UV 
convergent, IR convergent, and is completely indepen- 
dent of g, which can be scaled out by defining new vari- 
ables p' = p/q and e' = Imejq^. It follows that the the 
diagrams representing repeated scattering in the particle- 
particle and particle-hole channels do indeed cancel, so 
that t aaaaZ^ does not renormalize. 

Combining this with our argument demonstrating that 
abbaZ"^ does not renormalize at log^ order (see dis- 
cussion below Eq.(32)), and recalling that ^ aaaa — 
"^AABB, we conclude that we can set T = with log^ 
accuracy. 



VII. SOLUTION OF RG FLOW EQUATIONS. 
ZERO BIAS ANOMALY IN BILAYER 
GRAPHENE 

Since the only quantities which renormalize at log^ or- 
der in a one loop RG are the quasiparticle residue Z and 
the interaction vertex function F, the problem of find- 
ing the RG flow of these quantities reduces to solving 
Eas. (P71) . (P^ . All other quantities do not renormalize at 
log square order, and may thus be treated as constants 
with logarithmic accuracy. 

Additional simplification arises due to the Ward iden- 
tity TZ = 1, Eg. ([50)1 ■ Using it to decouple the RG equa- 



dZ 



a + c) Z, 



(44) 



where we retained a constant c = In ^° correspond- 
ing to the first term in the self energy renormalization, 
Eq.dSl. 

Integrating the RG equation, and taking into account 
the boundary conditions Z{Q) = F(0) = 1, we obtain 

^(O-xp(-^).F-(e), e^m^. (45) 

We note that in the limit of small £,'^/N, we reproduce the 
perturbative result^ for the residue, Ea. ([M)) . However, 
our result (^5]) applies for all ^, both small and large. The 
fermion propagator at arbitrary energies and momenta is 
then given by 



G{u,k) = -Z{0 



lUJ - 



i^o(k) 



(46) 



At zero temperature, the infrared cutoff is supplied by the 
external frequency and momentum, such that = In 
and A = + (A;2/2m)2 

Thus, the quasiparticle residue in undoped BLG is sup- 
pressed to zero by electron-electron interactions, Ea. (l46t . 
This is reminiscent of the situation in disordered metals, 
where enhancement of interactions by disorder produces 
a renormalization of electron self energy of a log^ formii, 
and analysis of an RG flow-^^ yields a suppression of the 
quasiparticle residue similar in form to our Eq.(|46p. The 
suppression of the quasiparticle spectral weight at low 
energies, governed by the Z{S^) dependence, will manifest 
itself directly in the behavior of the tunneling density of 
states of BLG, similar to disordered metals. 

We note parenthetically that, while keeping the con- 
stant term c in the RG equation (|44l) is formally beyond 
the log^ accuracy generally adopted in our analysis, it 
can be justified on the same grounds as in the discus- 
sion of the zero bias anomaly in disordered metalsi^ii^. 
Because of its fairly large value for = 4, given by 
c = In27r2 ss 2.98, this term may significantly alter pre- 
dictions for the behavior of Z at intermediate energies 

To analyze the suppression of tunneling density of 
states (TDOS), we use its relation to the retarded Green 
function^'*^ 



p[uj) = --Im [TrGfl,(w,k)] 

TT 



(47) 



where Gi^(cLi,k) is obtained from the Matsubara Green 
function analyzed above, Ea.(pS|. by the analytic con- 
tinuation of frequency from imaginary to real values, 
joj — > w + irj. 

It is convenient to take the trace before performing 
the analytic continuation. The trace may be most easily 
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taken in a basis of free particle eigenstates (plane waves 
with appropriate spinor structure), which amounts to in- 
tegrating Eq. (|46l) over all k values, Tr G = / G{lu, k)d'^k. 
Noting that the term containing i?o(k) vanishes upon in- 
tegration due to the angular dependence, we write 



TrG 



2Nn 



z{0- 



-dz. 



(48) 



where z — k^/2m and A^o is the density of electronic 
states in BLG in the absence of interactions. 

It can be seen that the integral over z is determined by 
z ^ uj. It is therefore convenient to introduce a variable 

if = sinh^"'^(z/cLi) and write 



TrG 



2No 



Z{S,u — In cosh (/j) 



dip 
cosh (fi ' 



(49) 



where = ln(Ao/ti;). Noting that this integral is dom- 
inated by (/3 ^ 1, we obtain an estimate of the spectral 
weight: 



(50) 



The form of this expression remains unchanged, to lead- 
ing log^ order, upon analytic continuation to real fre- 
quencies. 

The expression in Ea. ((5U)) can be re-arranged by using 
Eq.dini) as 



p{uj) = Nq exp - 



ln2 iv!^ _ 1^,2 N 



8Ao 



Ntt^ 



(51) 



Thus, we see that the only effect of the UV cutoff Aq is 
to rescale the prefactor for the TDOS without affecting 
the frequency dependence. Absorbing the dependence on 
Aq in the prefactor, we have 



p(w) = No exp - 



1 



In 



(52) 



Tunneling measurements yield p{uj = eV), where V is 
the bias voltage. The interaction suppression of the 
TDOS, Eq . (f50| . will therefore manifest itself as a zero 
bias anomaly in tunneling experiments. The predicted 
behavior the TDOS is shown in FiglJ] Because of the ex- 
ponential dependence in Eg. ((52)) . the suppression rapidly 
becomes more pronounced at lower energies. 

Closing our discussion of the zero bias anomaly in 
BLG, we note that the results described above apply only 
to the system at charge neutrality. Away from neutrality, 
with the Fermi surface size becoming finite, the effects of 
screening will grow stronger, resulting in a weaker ef- 
fective interaction. Yet, even in this case, the tunneling 
density of states will be described by the suppression fac- 
tor p{ijj = eV)/NQ given by Eg. ([SO)), provided that the 
bias voltage eV exceeds the Fermi energy measured from 
the neutrality point. 
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FIG. 3: Tunneling density of states (TDOS) of BLG at charge 
neutrality, Ea. (|52p . is shown as a function of external bias 
oj = eV. Predicted TDOS is shown for two different values 
of the dielectric constant in Eq, Eq.Q: k = 1 (solid curve) 
and K — 2.5 (dashed curve), describing free-standing BLG 
and BLG on SiO substrate, respectively. Plot is normalized 
so that p = 1 at an external bias of 100 meV. 



VIII. SINGLE LOG RENORMALIZATION OF 
ELECTRON MASS 

Thus far we have concentrated on log^ flows. However, 
the analysis may be extended to obtain the subleading 
single log flows of the action. We illustrate this procedure 
by calculating the renormalization of the mass (which did 
not renormalize at log^ order in the RG). This calcula- 
tion is interesting because it allows us to investigate the 
interaction renormalization of the compressibility-a di- 
rectly measurable quantity, and also because it allows us 
to illustrate how much slower the single log flows are than 
the log^ flows. 

In this section, we first analyze mass renormalization 
by extracting it directly from the self energy. After 
that, in Scc lIXI we consider electron compressibility of 
BLG and show that the log divergent correction to the 
compressibility matches exactly our prediction for mass 
renormalization obtained from the self energy. 

In BLG, the self energy is a 2 x 2 matrix, given by 
Eq.([T4|)), which is related to the renormalized Green 
function by the Dyson equation, 

G-l(c.,q) = Gol(c.,q) - ( ^^^J'^'q) ^ABiu^,n) 

(53) 

As discussed in Sec lIVl the leading log contribu- 
tion to the self energy is proportional to Gq^ , since 
dYjAB / d{q\/2m) = idYjAA/duj. This means that all 
renormalization can be attributed to the residue Z with 
mass remaining unchanged. However, as we now show, 
this equality is only true to leading logarithmic order. 

Comparison of Eq.lUS]) with Eq.(I15l) and Eq.([n]) indi- 
cates that the mass renormalization is given by 



5m 



m 



Zo 



AA 



AB 



did d{q\l2m) 



(54) 



Here, idYiAA/duj is defined by Eg. p^ . For the second 
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term, we obtain the expression 



51] 



AB 



d{q\/2m) 



ded p 
"(2^ 



^ 2m/ 



, 2m) 



^\2m) 



r^ZUie,p), 



(55) 



where U{e,p) is given by Eq.®. To evaluate the dif- 
ference in Ea. ((54)) . it is convenient to subtract the inte- 
grands of Eqs.( (|19l55p ) before doing the integrals. Once 
again, we use the "polar" representation of the frequency 
and momentum variables, lu — rcos9, p^/2m = rsin^, 
and obtain 

dm 
m 



dr /"^ de rgZ^ (3 sin^ 6* - 4 sin"* 
27r V2rsin6l- ^U{e) 



where n(^) is the polarization function introduced in 
Eq.(21), and r is measured in units of Eq as before. The 
integral over 9 is now fully convergent, and the resulting 
expression is only single log divergent. Integrating ana- 
lytically over r and then integrating numerically over 9, 
we find 



TO 2A^7rln4 ° ° A' 



(56) 



Converting this recursion relation into a differential equa- 
tion, we obtain 



din TO 

~dr 



0.56 



2iV7rln4 



(57) 



This equation cannot be solved for general ^ by apply- 
ing the Ward identity Eq. ipD]) . since the Ward identity 
only holds at leading log^ order, while the mass flows at 
subleading (single log) order in Eq. ((57|) . In the pertur- 
bative limit -^^ <C 1, when Z k, \ and F w 1, we obtain 
a logarithmic correction to the mass 



m(^) = to(0) 1 + 



0.56 



2iV7rln4 



(58) 



We may relate this mass renormalization to a measur- 
able quantity, by noting that the electronic compressibil- 
ity K is proportional to the density of states which is 
proportional to the mass. Thus, the logarithmic renor- 
malization of the mass in Eq. (|58|) should manifest itself 
in a logarithmic enhancement of the electronic compress- 
ibility. The relation between mass renormalization and 
compressibility will be further discussed in Sec lIXI 



IX. INTERACTION CORRECTION TO 
COMPRESSIBILITY 

Here we explicitly calculate the renormalization of the 
compressibility. By doing this we shall confirm that 



the compressibility does not renormalize at leading (log 
square) order, and also extract the single log renormal- 
ization of the compressibility. The interaction correction 
to the compressibility K is given by 

where /i is the chemical potential, and F is the interac- 
tion energy. Within the RPA framework, the interaction 
energy is expressed as 



dujd^p 
"(2^ 



ln(l-y(q)n(Ai,^,q)) 



(60) 



Here, Il{fi,uj,q) is the non-interacting polarization func- 
tion evaluated at a chemical potential /i, and V{q) is the 
unscreened Coulomb interaction V{q) — l-Ke^ j Kq. 

To evaluate the second derivative in (15^1) . we consider 
the difference AF = F(/i) — F(0). After rearranging logs 
under the integral, we rewrite this expression as 



Ai^ = - 



(27r)3 



In l-C/^,,(H(^,c^,(Z)-H(0,c<.,q)) 



(61) 

where now Ui^,q is the dynamically screened Coulomb in- 
teraction, Eq.(l5|). Since the compressibility is obtained 
from the free energy through K — —d^F/dfj?, the prob- 
lem of calculating the interaction renormalization of the 
compressibility is reduced to that of calculating the po- 
larization function at finite /x. This may be calculated 
through methods similar to those developed in Refl^. We 
define e± = e ± cij/2, p± = p ± q/2 and z± — |p±p/2m. 
The non-interacting polarization function at finite /i is 
given by 

n(/i,w,(j) = TrG'(Ai,e+,p+)G(Ai,£-,P-) 
ded^p 1 



Tr 



2N 



(2^)3 {ie+ - /i - Hoip+)) {is- - /i - Ho(P-)) 
ded'^p 

(27r)3 {e+ +i{ji + z+)) (e+ + i{n - z+)) 
{ie+ — f^){is- — fi) + z^z^ cos 29pq 
(e_ + i(/i + z_)) (e_ + i{fi - z_)) 



(62) 



where 9pq is the angle between p_|_ and p_. We now 
perform the integral over e by residues to obtain 



H(/i,w,q) =iV 



(Pp {z+ 



cos 29pq)Q{z+ — fi) 



N 



+ (w,q- 
'+=^ d^p^ 



{2nY 



■ oj^ -|- 2iwz+ 



(27r) 



-q) 

1 



2z_ 



Jpq 



-f (w, q 



-w,-q) 



(z+ + iujY — z_ 



(63) 



In the limit /i — ^ 0, this reproduces the non-interacting 
polarization function from Ref.Q- Now we expand 
Eg. dFlj) to leading order in small ^ to obtain 



AF 



^U{u^q)^^. 



(64) 
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The term linear in /i must vanish, by particle hole sym- 
metry. Taking derivatives of Eg. ([55)1 greatly simplifies 
the calculations, since it turns the two dimensional inte- 
gral over momenta into a one dimensional integral over 
momentum angles, which is fully convergent, and may be 
evaluated numerically. We find 



(w2 + ^2)2 ' -1 

^[/(w,q) 



2to' 



(27r)3 



^^J? 



(65) 
(66) 



We again change to the coordinates lu — rcosO, Zq = 
rsmO, and measure r in units of Eq. Note that even 
though the interaction has a pole at 6* — > 0, vr, this pole 
is canceled by d^Ii/dfi^ having a zero at 6* — 7> 0, tt. As 
a result, the 9 integral is fully convergent. Integrating 
numerically over 9 and analytically over r, we find that 
the fractional change in the compressibility is 



KiO) 



0.56 

2A^7rln4 



(67) 



a result that agrees exactly with Eq. ([58|) . We note that 
an enhancement of the compressibility due to interac- 
tions was also predicted in Ref.|19j. However, the effect 
described by Eg. ([55)) is much weaker than that predicted 
in Ref.fl^, because we have worked with a screened in- 
teraction, whereas in Ref.[l9j screening was not taken 
into account. 

In summary, the compressibility does not renormalize 
at leading (log square) order, just as in the Luttinger 
liquids, and while there is a subleading logarithmic cor- 
rection, the pre-factor is quite small (0.56/(2Af7rln4) « 
0.016). Thus, in contrast to the zero-bias anomaly in 
TDOS, experimental detection of the interaction correc- 
tion to the compressibility is likely to be challenging. The 
difference arises because the single log flows are much 
weaker than the log^ flows, retrospectively justifying our 
earlier neglect of the single log flows in the RG. Hence, 
strong suppression of the tunneling density of states at 
energy scales where the compressibility is not signifl- 
cantly renormalized is a key signature of the marginal 
Fermi liquid physics in bilayer graphene. 



X. DISCUSSION AND CONCLUSIONS 

Here we briefly discuss the range of validity of our re- 
sults. Our analysis was organized as a perturbation the- 
ory in r^Z^/TV. Since TZ = 1 at leading (log square) 
order, the perturbation theory remains well defined un- 
der the log square flows. However, our analysis neglected 
subleading single log flows. For ^ « TVtt^, the sublead- 
ing single log flows become important, and the analysis 
leading to the expression Eg. ps)) no longer applies. A 



mean fleld theory of subleading single log effect^ indi- 
cates that a gapped state develops at ^ = ^A^tt^, the 
scale which we tentatively identify as the limit of validity 
of our analysis. 

How can the marginal Fermi liquid physics be distin- 
guished from the formation of a gapped state? We note 
that at very low energies, once the gapped state has de- 
veloped, the tunneling density of states will vanish any- 
way. However, in the gapped regime, the compressibil- 
ity will vanish also. What we have shown, however, is 
that there is a large range of energies greater than the 
energy scale for gap formation, where the tunneling den- 
sity of states vanishes, while the compressibility remains 
essentially unchanged. Such behavior represents the key 
signature of the marginal Fermi liquid physics discussed 
above, which is analogous to the Luttinger liquid physics. 

In our analysis, we neglected the short range inter- 
actions which are characterized by lattice scale, such 
as the interlayer density difference interaction V- — 
\{Vaa — Vab) — Tre^d and the Hubbard-type on-site re- 
pulsion. Short range interactions are non-dispersive, do 
not renormalize the Green function in the weak coupling 
limit, and hence do not alter our results. Short range in- 
teractions also produce only single log renormalization^ 
and therefore do not need to be included in our log square 
RG. Similarly, we justify our neglect of the trigonal warp- 
ing effect^ by noting that trigonal warping is significant 
only on energy scales smaller than the characteristic en- 
ergy scale for onset of gapped statesS'. 

Finally, we note that our analysis made use of the fact 
that there were no un-canceled log square divergences at 
one loop order in the RG, except for the renormaliza- 
tion of the quasiparticle residue and the Coulomb vertex 
function, which were related by a Ward identity, Eq. (l30t . 
Technically, in order for our neglect of higher loop cor- 
rections to be justified, we also require that there are no 
un-canceled log square divergences beyond one loop or- 
der in the RG, except those that are constrained by Ward 
identities. We believe this to be the case, however, the 
proof requires a non-perturbative approach, which lies 
beyond the scope of the present work. 

To conclude, we have examined the one-loop RG flow 
for bilayer graphene. We have demonstrated that the 
quasiparticle residue Z and the Coulomb vertex function 
F both flow as where ^ is the RG time. All other 
quantities flow only as ^. The structure of the RG for 
Coulomb interacting BLG has strong similarities to the 
RG for the one dimensional Luttinger liquids. In partic- 
ular, we predict a strong interaction suppression of the 
tunneling density of states for undoped BLG, even at en- 
ergy scales where the electronic compressibility is essen- 
tially unchanged from its non-interacting value. These 
predictions may be readily tested by experiments. 

We acknowledge useful conversations with A. Potter 
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